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Gaussian quantum states hold special importance in the continuous variable (CV) regime. In 
quantum information science, the understanding and characterization of central resources such as 
entanglement may strongly rely on the knowledge of the Gaussian or non-Gaussian character of the 
quantum state. However, the quantum measurement associated with the spectral photocurrent of 
light modes consists of a mixture of quadrature observables. Within the framework of two recent 
papers [Phys. Rev. A 88, 052113 (2013) and Phys. Rev. Lett. Ill, 200402 (2013)], we address 
here how the statistics of the spectral photocurrent relates to the character of the Wigner function 
describing those modes. We show that a Gaussian state can be misidentified as non-Gaussian and 
vice-versa, a conclusion that forces the adoption of tacit a priori assumptions to perform quantum 
state reconstruction. We experimentally analyze the light beams generated by the optical parametric 
oscillator (OPO) operating above threshold to show that the data strongly supports the generation 
of Gaussian states of the held, validating the use of necessary and sufficient criteria to characterize 
entanglement in this system. 
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I. INTRODUCTION 


Entanglement and non-Gaussian quantum states are 
prized resources in the field of quantum information with 
continuous variables (CV). While Gaussian states per se 
allow the generation of highly entangled states and thus 
the realization of important quantum information proto¬ 
cols 0,1, certain tasks such as quantum computation 
and entanglement distillation require non-Gaussian fea¬ 
tures, and thus the ability to correctly recognize the avail¬ 
ability of those resources in actual physical systems 11- 

Light is perhaps the most important physical system 
for CV quantum information. Measurement of the quan¬ 
tum state of light is realized by the statistical analysis 
of photocurrent fluctuations, the quantum noise [51-41 lj]. 
The physical objects of interest are often field modes with 
well defined frequency, in which case the spectral noise 
power of photocurrent fluctuations is used to provide in¬ 
formation about the field quantum state. 

However, owing to a lack of phase coherence be¬ 
tween the quantum state and the spectral components 
of the photocurrent signal, current measurement tech¬ 
niques only provide a pure quantum measurement for 
a restricted class of quantum states possessing a strong 
degree of symmetry [12l4l4| . In most experiments, the 
quantum observable associated with the spectral pho¬ 
tocurrent provides a mixed measurement of the two side¬ 
band modes and requires the adoption of a priori as¬ 
sumptions to interpret the data in terms of moments of 
field quadrature operators. These limitations constitute 
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a problem for the demonstration of quantum informa¬ 
tion protocols with truly general and unknown quantum 
states of spectral light modes. 

Yet, it is generally believed that Gaussian spectral pho¬ 
tocurrent statistics can be taken as proof of the Gaus¬ 
sian character of the quantum state. Conversely, the ob¬ 
servation of non-Gaussian photocurrent statistics is ac¬ 
cepted as strong evidence of non-Gaussian features of the 
field quantum state [lfj. Contrarily to this belief, we 
show that the mixed character of the photocurrent mea¬ 
surement operator fundamentally masks the phase space 
statistics of field quadratures, implying that a Gaussian 
state may appear non-Gaussian and vice-versa. Such du¬ 
biety about the statistics obeyed by the CV quantum 
state is especially harmful to experimental quantum in¬ 
formation, where non-Gaussian properties of field modes 
stand as a necessary resource to perform more powerful 
tasks: misidentifying those resources leads to incorrect 
or, at least, unreliable implementations. 

In this paper, we perform a formal analysis of the mea¬ 
surement operator moments associated with the spec¬ 
tral photocurrent to evaluate the extent to which the 
observed statistics can be interpreted in terms of quan¬ 
tum state features. We consider in Sec. eh the photocur¬ 
rent fluctuations as an incoherent mixture of two inde¬ 
pendent quantum measurements to show, by exploiting 
higher-order moments, that a Gaussian photocurrent can 
be produced either by Gaussian states with spectral two¬ 
mode symmetry or by very peculiar (and hence implausi¬ 
ble in most experiments) non-Gaussian quantum states. 
In fact, our reasoning links the Gaussian character and 
the symmetry of the two-mode field as equivalent a pri¬ 
ori knowledge always needed in the usual experimental 
setup (but rarely mentioned) to reconstruct the spectral 
quantum state. We additionally point out that the obser- 
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vation of non-Gaussian spectral photocurrent can not be 
readily associated with a non-Gaussian quantum state, 
requiring further investigation and possibly the improve¬ 
ment of the quantum measurement technique of CV spec¬ 
tral modes to a new level of experimental rigor. 

We then employ our methods in Sec. m to experi¬ 
mentally investigate the quantum state produced by the 
optical parametric oscillator (OPO) operating above the 
oscillation threshold, establishing with great confidence 
that the associated photocurrent is indeed Gaussian and 
indicates a Gaussian quantum state displaying modal 
symmetries in the spectral modes linked by the paramet¬ 
ric process. This conclusion substantiates the use of the 
semi-classical treatment in the above-threshold OPO to 
effectively halve the number of relevant modes and apply 
necessary and sufficient criteria for multipartite entan¬ 
glement directly to the spectral matrix [f6l-[l8|. It is our 
feeling that such an experimental demonstration fills a 
hole that went unrecognized for some time. 

Finally, to guarantee that the measurement operator of 
the spectral photocurrent yields the statistics of pure field 
quadrature observables, we offer our concluding remarks 
in Sec. EY1 where the use of phase coherent optical and 
electronic local electronic oscillators is proposed to realize 
a doubly phase-coherent detection. Only in this improved 
scenario can the quantum noise achieve the status of a 
formal pure measurement operator in the CV regime of 
spectral modes of light. 

II. GAUSSIAN PHOTOCURRENT AND 
PROPERTIES OF THE SPECTRAL FIELD 
MODES QUANTUM STATE 

Spectral photocurrent operator - The quantum state 
of spectral (sideband) modes is measured by the pho¬ 
tocurrent observable [191 ], described by the measurement 
operator /( t ) and obtained by the beating of the optical 
local oscillator (LO) and the field modes of interest po| . 
Any Fourier component of this signal can be extracted 
by mixing it with a sinusoidal electronic reference - the 
electronic local oscillator (eLO) - at frequency Cl, yield¬ 
ing the spectral photocurrent operator Iq defined as 



The operator In(t) represents the (complex) amplitude 
of the photocurrent beatnote signal measured at time t. 
The integration limits are determined by the spectral res¬ 
olution A Cl = 2n/T of the electronic downmixing process 
and define the spectral width of measured upper ujq + Cl 
and lower uiq — Cl sideband modes. The spectral pho¬ 
tocurrent operator In(t) gathers in one quantity the two 
observables associated with the photocurrent cosine I cos 
and sine J s j n components, given by 

Iq = (/cos + */sin )• (2) 
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FIG. 1. Scheme to measure electronic quadrature components 
of each photocurrent signal. The photocurrent is mixed with 
two electronic references in quadrature. 

We may interpret I CO s(t) and / s i n (t) as two independent 
single quantum measurements taken at time t [14j . Fig. 
[T] illustrates the measurement apparatus. These observ¬ 
ables access the quadrature operators of optical sideband 
modes at frequencies ujq ± Cl, where ujq is the LO optical 
frequency, and bandwidth A Cl. These longitudinal modes 
involve the annihilation operators a± n = (p±n + iq±n)/ 2, 
where p±o, and q±n are respectively the field amplitude 
and phase quadrature observables satisfying the canon¬ 
ical commutation relation [pn,qn'] = 2 iS(Cl — Cl'). The 
interpretation of the photocurrent observables I cos and 
/ s i n in terms of the spectral quantum state is easily de¬ 
scribed after a change of modal basis by the symmetric 
(5) and anti-symmetric ( A ) combinations of sideband 
modes, defined as 

hs — 77^(ttajQ-i-o T h Uo — q), 

ha, — "7/2 (^<-^0+^ hug — o). ( 3 ) 

The exact connection between the quadrature observ¬ 
ables of these modes, defined by the expressions a s = 
p s + iq s and a a = p a + iq a , and the photocurrent observ¬ 
ables I cos and / s in depends on the measurement tech¬ 
nique in use. In homodyne detection (HD), the cosine 
observable is a pure quantum measurement of the form 
I C os = -/cos (Ps, q s ), where L cos represents a linear com¬ 
bination; the sine observable has the same general form, 
but depends only on quadratures of the A mode. In res¬ 
onator detection (RD), both observables combine S and 
A modes in independent ways, so that / cos and I s , n carry 
information on both modes [l4j . 

Phase mixing and the quantum state - Measurement 
of the photocurrent components / cos and / s i n requires a 
well defined phase relation between LO and eLO, a condi¬ 
tion in general not satisfied in experiments. In fact, that 
limitation does not affect each individual quantum mea¬ 
surement of the spectral photocurrent, typically carried 
out on a time scale much shorter than the characteristic 
time of LO phase diffusion; they realize the pure quadra¬ 
ture observable 

Ig = cos 8 Icos + sin 8 J sin . (4) 

However, the collection of quantum statistics requires 
many such single measurements and hence a time inter- 
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val that greatly surpasses the relative coherence time be¬ 
tween LO and eLO. As a result, photocurrent moments 
involve averaging different quadrature directions in phase 
space, entailing mixed quantum statistics (^-averages of 
moments of Ig on the quantum state). Let us construct 
the operator Slg = Ig — (Ig) with zero average on the 
quantum state, and consider its 0-average variance, de¬ 
noted as A 2 Ig. Using Eq. (0}, 

A 2 ig = 3 hf d0'((6ig + g,) 2 ) 

= ±A 2 / cos +±A 2 / sin , V0, (5) 

a quantity actually independent of the choice of 0 (i.e. 
no LO-eLO relative phase information remains, as ex¬ 
pected). Furthermore, components in quadrature (SIg 
and SIg + tt/ 2 ) are uncorrelated, 

(Sig5Ig + „/ 2 ) = ^ J dd'(Sig + g'SIg +n / 2 +g') 

= (SIcos SIsin - SI sin SI cos ) = 0, V0, (6) 

since [SIg, SI e+ ^/ 2 ] = 0. In this so-called phase mixing 
regime , in which most experiments dealing with the spec¬ 
tral quantum noise of light are realized, Eqs. © and © 
show that any electronic component of the spectral pho¬ 
tocurrent displays the same quantum statistics and is un¬ 
correlated with its simultaneously measured (in quadra¬ 
ture) component: in other words, the measured spectral 
photocurrent is stationary [u|. But since stationarity 
holds simply as a consequence of phase mixing, it does 
not convey any information on the quantum state; in 
fact, all such information lies in the variance A 2 Ig, the 
quantity proportional to the spectral noise power S(Cl) 
of semi-classical models, 

S(Sl) = (SinSi-n) = |A 2 / cos + ±A 2 / sin . (7) 

That is probably the reason why in most experiments 
the spectral photocurrent is measured directly on a spec¬ 
trum analyzer without further concern about the pure 
measurement operators I cos and I s i n , give n that their 
statistical properties seem inaccessible 21]. Neverthe¬ 
less, the current approach to the measurement of spectral 
moments as solely based on the mixed moment of Eq. © 
can not be considered a pure quantum measurement for 
most quantum states [lj, H3|. Furthermore, it fails to 
provide all the information required to reconstruct the 
two-mode spectral quantum state. Let us consider Gaus¬ 
sian quantum states (i.e. the simplest class of interesting 
quantum states) to illustrate this point. The complete 
reconstruction of any quantum state in this set would 
require the measurement of 10 independent second-order 
moments, gathered in the two-mode covariance matrix as 

/ A 2 p s C(p s q s ) C(p s p a ) C(p s q a ) \ 

_ | C(p s q s ) A 2 q s C(p a q s ) C(q s q a ) | , , 

C(p s Pa ) C(p a q s ) A 2 p a C(p a q a ) ’ ^ 

V c{p s q a ) C(q s q a ) C(p a q a ) A 2 q a ) 


where quadrature operator variances are denoted as in 
e.g. A 2 p s = ((p s — (ps)) 2 ) and symmetrized correlations 
as in e.g. C(p s q s ) = ((p s -(p s ))(q s -(q s )) + (q s -(q s ))(ps- 
(p s )))/2. It is clear that measurements of S( fl) defined 
in Eq. © can provide only a fraction those 10 moments, 
although the exact amount of information is dependent 
on the particular measurement technique in use. HD 
can only retrieve 3 combinations of moments, while RD 
accesses 1 additional combination [lij j. However, if one 
could see beyond the phase mixing process to establish 
the phase-locked photocurrent components I cos and I sln 
as stationary, 6 combinations of moments could be quan¬ 
tified as null and the covariance matrix would simplify 
to 


V = 
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(9) 


In this case, HD would be able to access the local 
single-mode sector of the quantum state (either S or 
A), and RD would provide complete two-mode quantum 
state reconstruction (by additionally accessing the S mo¬ 
ment related to the energy asymmetry of spectral modes 

±o) g. 

Hence establishing the stationarity of the phase-locked 
photocurrent components I cos and I s - ln (or tacitly assum¬ 
ing it) stands as a central part of spectral quantum state 
measurement, either in the single-mode picture (HD) or 
in the complete two-mode picture (RD). We show next 
how it is possible to ascertain with great confidence the 
stationary property even in the presence of phase mixing, 
by analyzing higher order moments of the spectral pho¬ 
tocurrent. The reasoning we present provides the neces¬ 
sary conceptual basis to correctly interpret the photocur¬ 
rent quantum noise in terms of moments of the symmetric 
covariance matrix of Eq. ©, clarifying the tacit assump¬ 
tions employed in nearly all experiments with spectral 
modes of light in the CV regime. 

Gaussian character of the measured photocurrent - We 
investigate next what can be said about the two-mode 
quantum state by having access only to the mixed oper¬ 
ator moment of Eq. ©• Particularly, we ask what can be 
learned about the quantum state when the photocurrent 
is faithfully established by the experiment as Gaussian or 
non-Gaussian. Let us consider how higher order moments 
of the photocurrent relate to higher order quadrature op¬ 
erator moments in the phase mixing regime. 

We denote the (phase mixed) 2?r th -order photocurrent 
moment as cr^ 2n \ n = 1, 2,3,..., defined as 

* {2n} = h r m) 2n )do. (io) 

J —TT 

These measured moments can be expressed in terms 
of quadrature operator moments (and therefore related 
to the two-mode quantum state), through the measure¬ 
ment operators Sl cos and Sl s ; n of Eq. ©. These observ¬ 
ables correspond to what would be measured with doubly 
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phase-coherent detection, and represent the phase-locked 
photocurrent components masked by the phase mixing 
as well as pure independent measurements of the spec¬ 
tral quantum state. In order to relate those quantities 
to the mixed photocurrent moments of Eq. m, we de¬ 
note the moments of the phase-locked photocurrent as 
4os" } = ((d/ C os) 2 "} and cr s { ; 2 " } = (( SI sin ) 2n }. 

For the particular but important case of Gaussian 
statistics (applicable to both the quantum state and the 
measured photocurrent), odd moments (other than the 
first) are null and all even moments can be expressed in 
terms of the variance of the distribution. This fact mo¬ 
tivates us to denote the stand ard deviations of distribu¬ 
tions as follows: s = y/ crl 2 1 for the measured (mixed) 

photocurrent, and s cos = \jcrial or s sin = \J o4/ 
for the cosine or sine components subjacent to phase 
mixing. We also denote their normalized correlation 
(equal to the correlation between spectral modes) as 
C = (d/ C os^d s i n )/(s CO sS s i n ) ( 1 ^ C ^ 1). 

We would like to establish how the Gaussian character 
of the photocurrent constrains the field quantum state. 
Let us consider for example the fourth-order moment 
cj{ 4 } of the mixed photocurrent. On the one hand, a 
Gaussian photocurrent would imply cr^ = 3s 4 ; on the 
other, Gaussian quantum states would require fourth- 
order moments fulfilling = 3s 4 os , er^ = 3s 4 in and 

(Wcos^sin) = (! + 2c 2 )s 2 os s 2 in . 

To better analyze the actual distributions in terms 
of the Gaussian statistics, we define ‘deviations’, i.e. 
quantities constructed to be null for Gaussian distri¬ 
butions, as follows: 6 = — 3s 4 for the measured 

photocurrent and, for the quantum state, the quanti¬ 
ties <5 cos = cJos - 3s 4 os , d sin = c>4 } - 3s 4 in and S c = 

{Sic :os<4L) - (! + 2c 2 )Scos4n- 

We then employ Eqs. © and ed to relate the (fourth- 
order) photocurrent deviations and quantum state devi¬ 
ations from the Gaussian statistics by 

= 5 cos + dsin + 2 S c + (sL — S 2 in ) + 4c 2 S 2 os S 2 in . (11) 

Interpretation of Eq. ED requires careful consideration 
of possible experimental scenarios. We note that actual 
measurements can only access the left side of the equal¬ 
ity above with the goal of determining all the terms on 
the right: an impossible task if no a priori assumptions 
are allowed. Let us analyze the exact content of those 
assumptions in the usual spectral noise detection. 

The first possibility is that the experimental data es¬ 
tablishes the photocurrent statistics as Gaussian, impos¬ 
ing the constraint <5 = 0, i.e. 

^cos + 5 s in + 25 c = — (s cos — S s in) — 4c S cos S s ; n . (12) 

This equation establishes that infinitely many quantum 
states can in principle give rise to Gaussian spectral pho- 
tocurrent. However, their fourth- and second-order mo¬ 
ments must fulfill a very stringent relation. Hence most 


quantum states satisfying Eq. ED- even though not ruled 
out by the data alone, would be very unusual, obliging us 
to consider different scenarios to help interpret the data. 

In the first scenario, we note that the quantum dynam¬ 
ics capable of producing such stringent quantum states is 
unlikely to be taking place in most experiments in quan¬ 
tum optics or atomic physics. These type of experiments 
usually employ at least one bright light beam to par¬ 
ticipate in the quantum dynamics, in which cases a lin¬ 
earized or mean field approach in general supports the 
onset of Gaussian quantum states. The observation of 
Gaussian mixed photocurrent [Eq. (fl2l) ] together with 
the assumption of Gaussian quantum states, implies the 
new condition 

(«cos - «sin) 2 + 4c 2 S 2 os 4n = 0, (13) 

which can only be fulfilled if 

Scos — ^sin and c — 0, (14) 

meaning that the cosine and sine photocurrent compo¬ 
nents must present the same variance and be uncorre¬ 
lated: a statement of stationarity [l4j . Even though 
phase mixing always leads to stationary mixed photocur¬ 
rent, verifying that its fourth-order moment is compat¬ 
ible with the Gaussian statistics makes it a better case 
to establish as stationary the phase-locked photocurrent 
related to pure quantum measurements. Thus checking 
higher order photocurrent moments allows us to remove 
some hindrances of the phase mixing process and par¬ 
tially ‘see through it’. 

As noted previously, the most relevant consequence of 
stationarity in our context can be established by inter¬ 
preting Eq. ED in terms of properties of the spectral 
quantum state. Assuming it as Gaussian, stationarity 
implies the form of Eq. © for the two-mode quantum 
state, which in turn allows us to employ only the mixed 
photocurrent variance as source of information about 
the quantum state. Only after this formal verification 
can one confirm that both HD and RD are able to ac¬ 
cess the single-mode quantum state of either mode S or 
A, justifying the application of entanglement criteria di¬ 
rectly to the spectral noise matrix, equal in this case to 
the single-mode covariance matrix of the quantum state 
(semi-classical approach). If necessary, RD can be used 
to further access exclusive two-mode spectral features, a 
capability especially important in case one needs to re¬ 
construct not only the local quantum state of S or A 
mode, but also the local quantum states of individual 
sideband modes ±H, in which case the energy imbalance 
becomes essential to perform the change of modal basis 
on the two-mode quantum state. 

We may summarize the first scenario that may be 
adopted to interpret the observation of Gaussian mixed 
photocurrent as employing the following set of concepts 
in a self-consistent manner: Gaussian quantum state -o- 
Stationary phase-locked photocurrent -o- Symmetric S 
and A local quantum states. If one of these concepts 
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is adopted as a priori assumption to help interpret the 
experimental data (since none of them can be directly 
checked due to measurement mixedness) the others fol¬ 
low as consequences. For instance, in the case of phe¬ 
nomena associated with parametric downconversion, the 
broadband nature of the physical interaction will proba¬ 
bly favor the adoption of ‘symmetric S and A local quan¬ 
tum states’ as assumption, from which follow as conse¬ 
quences the validity of the semi-classical approach and 
the fact that an EPR-type quantum state is produced on 
the spectral sidebands (by a change of modal basis). 

In the second scenario, Eq. m allows the possibility 
that Gaussian photocurrent could be produced by non- 
Gaussian quantum states. Additional evidence would 
certainly be required in this case, since a Gaussian pho¬ 
tocurrent would hardly have any credibility in attesting 
the non-Gaussian character of the quantum state. The 
type of quantum state capable of fulfilling Eq. m would 
have to present very unusual connections between fourth- 
and second-order moments. As we show later, the con¬ 
ditions to be fulfilled become more stringent as moments 
of even higher degree are considered, making it a diffi¬ 
cult case for this kind of interpretation without further 
evidence. 

Since non-Gaussian features are desirable as a resource 
for quantum information protocols in CV, we consider 
the third scenario in which non-Gaussian mixed pho¬ 
tocurrent could be used as misleading evidence of non- 
Gaussian features of the quantum state. In fact, accord¬ 
ing to Eq. CEB). Gaussian quantum states impose upon 
the 4 th -order moment of the photocurrent the condition 

§<S = ( S cos - 4n) 2 + 4 C 2 4os4n- (!5) 

This identity states very generally that whenever the two 
individual S or A single-mode quadratures follow differ¬ 
ent Gaussian distributions, i.e. (( 8I cos ) 2 } ¥= ((<5/ s in) 2 ), 
the photocurrent will present non-Gaussian statistics. 
The photocurrent deviation from Gaussian statistics 
could thus follow as a consequence of mixed Gaussian 
processes, in fact the favored explanation in view of the 
ubiquity of Gaussian states in quantum optics. Such 
cases could be spotted by investigating the complete sta¬ 
tistical distribution of measured quantum fluctuations, 
wherein Gaussian mixed models could be employed to 
estimate the single-mode probability distributions subja¬ 
cent to the mixing of the measurement process. Failure 
to identify mixed Gaussian distributions would support 
further investigating the non-Gaussian character of the 
quantum state and probably require additional experi¬ 
mental capabilities. Conversely, Eq. m makes the ob¬ 
servation of Gaussian photocurrent as a special event, 
since its occurrence is a strong indication of Gaussian 
quantum states possessing modal quantum state symme¬ 
try. 

The conclusions reached for the 4 th -order moments are 
accentuated by extending our analysis to moments of 
even higher degrees. We find that the observation of 
Gaussian photocurrent in all orders either substantiates 


the Gaussian character of the quantum state or imposes 
increasingly stringent criteria upon the moments of the 
non-Gaussian quantum states capable of producing it. 
We define the deviation ^f 2n l of the 2n th -order moment 
from Gaussian statistics as 

j{2n} = a {2n} _ ^ _ jjj, ( 16 ) 


Similarly, the deviation <5{2n,2fc} from Gaussian statistics 
regarding the quantum state is defined as 


£ {2 n, 2fc} = 


(2 (n — k))\(2k)\ 2 (n-k) 2k 

2 n (n — k)\k\ cos sin ’ 
(17) 


where k = 0,1, 2,..., n and s = s cos = s s j n . Perform¬ 
ing the integral in 8 , Eq. m applied to higher-order 
moments yields 


) n 1 

{2n} _ n - ST' _^_ i/<rf \2ntcj \2(ra-fc)\ 

a cos) ( sm) h 

(18) 

In the following we suppose uncorrelated single modes 
(i.e. c = 0), a simplifying restriction that does not affect 
the physical conclusions. The connection between pho¬ 
tocurrent statistics and quantum state statistics provided 
by any even moment of order larger than 2 is 


j{2ra} 


n - 1 £ 

(2 n)\^{n-k)\k\ W- fc )’ 2fc > - 


(19) 


(2n- 1)!! 


k 


2 (n-k) 2k 
cos ^sin 


where the coefficients are d n .k = ,i! ( 2 (» 1 ) !! ^ 
The left-hand side of this expression relates only to pho¬ 
tocurrent deviations from the Gaussian statistics, and the 
right-hand side only to quantum state deviations. In fact, 
j{ 2 n} |-j ie on ]y parameter in Eq. (TTiil) available in the 
data. The quantum state must be determined from in¬ 
complete information due to the mixing process affecting 
the measurement operator. This is indeed an impossible 
task if any quantum state can be expected in the experi¬ 
mental dynamics under study. However, as noted before, 
usually that is not the case, and in fact the quantum 
dynamics producing non-Gaussian states constrained by 
Eq. m f° r a given measured value of df 2 "f is very diffi¬ 
cult to attain in the laboratory. Such state of affairs calls 
for a detailed analysis of Eq. m in particular cases of 
interest. 

In case the photocurrent is Gaussian, i.e. 8^ 2n ^ = 
0,Vn, Eq. (1191) provides very stringent restrictions con¬ 
necting the quantum state deviations from the Gaussian 
statistics with second-order operator moments. Those 
types of quantum states can usually be ruled out based 
on physical considerations of the quantum dynamics act¬ 
ing on the system. It is thus more plausible in this 
case to interpret the Gaussian photocurrent as indicat¬ 
ing the Gaussian character of the quantum state, for 
which the left-hand side of Eq. m must vanish (i.e. 
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^2(n—k),2k = 0,Vn, k). Assuming then the quantum state 
as Gaussian, one finds that the right-hand side of Eq. m 
can only vanish (as required to uphold the identity) if 
s C os = s s i n , in which case it is simple to verify that the 
sum of coefficients yields J2k^ n , k = 0- We find once 
more the self-consistent connections between Gaussian 
quantum states, symmetric A and S modes, and station¬ 
ary of the phase-locked photocurrents subjacent to the 
mixing process. Conversely, a Gaussian quantum state 
will produce non-Gaussian photocurrent with respect to 
every higher-order moment if s cos ^ s s j n . 

It is a striking fact that Gaussian photocurrent to¬ 
gether with the assumption of Gaussian quantum state 
will always lead to the equality of the single modes A 
and S quadrature distributions, and vice-versa. Similarly 
simple conclusions can not be reached for non-Gaussian 
photocurrents, since they could result either from truly 
non-Gaussian quantum states or from a mixture of two 
different Gaussian processes (i.e. when both single-mode 
quadrature distributions of mode S and A are Gaussian, 
but different). 


III. EXPERIMENT AND DATA ANALYSIS 

Experimental setup - The experiment (Figj2j) employs 
a frequency-doubled diode-pumped Nd:YAG cw laser 
(Innolight Diabolo) at 532 nm wavelength to pump a non¬ 
degenerate triply resonant optical parametric oscillator 
(OPO). Prior to injection in the OPO, the cw pump laser 
is filtered in transmission by an optical resonator (band¬ 
width of 0.7 MHz) to attenuate classical noise for analy¬ 
sis frequencies above 15 MHz. In this manner, the pump 
laser shows nearly vacuum noise limited quantum fluctu¬ 
ations at the spectral modes of interest at 12 = 21 MHz. 
The OPO consists of a type-II phase matched KTP 
(Potassium Titanyl Phosphate, KTi0P04) crystal with 
length l = 12 mm and a linear resonator with spheri¬ 
cal mirrors with 50 mm radius of curvature. The input 
coupler mirror has reflectivity of 70% for the pump (532 
nm) and > 99.9% for the signal and idler beams (~ 1064 
nm), and the output coupler mirror has transmission of 
96% for the infrared beams and is highly reflective for 
the pump beam. The OPO has a free spectral range of 
about 5 GHz and cavity finesses of 18, 135, and 115 for 
pump, signal and idler modes, respectively. The KTP 
crystal temperature is actively stabilized at 23°C. The 
OPO oscillation threshold power is Pth = 60(3) rnW. 
Operating above the threshold, it generates collinear sig¬ 
nal and idler beams with orthogonal polarizations, which 
are then spatially separated by a polarizing beam-splitter 
(PBS). The reflected pump beam is separated from the 
input beam by a circulator employing a Faraday rotator 
(RF) and a PBS. 

We employ the RD technique to reconstruct the quan¬ 
tum states of spectral modes in the pump, signal and idler 
beams individually [13j]. In addition to being a complete 
quantum measurement for Gaussian quantum states, in 


our experiment RD has also the advantage of allowing 
the use of the bright beams produced by the OPO as 
independent LO’s at the appropriate optical frequencies. 
One analysis resonator is employed to perform RD in 
each beam. Spatial mode matching between the beams 
and the analysis resonators is higher than 95%. The anal¬ 
ysis resonator labeled j = 0, corresponding to the pump 
mode, has spectral bandwidth of ss 12 MHz; resonators 
labeled j = 1 and j = 2, respectively used to measure 
signal and idler beams, have bandwidths of ~ 14 MHz. 
With these values it is possible to access all quadratures 
for the chosen sideband modes at the chosen analysis fre¬ 
quency of 21 MHz. 

Photodetection is performed with high quantum effi¬ 
ciency (> 95%) photodiodes (Epitaxx ETX300 for the 
twin beams and Hamamatsu S5973-02 for the pump 
beam). Quantum noise of each beam is measured with a 
pair of amplified photodetectors using the balanced de¬ 
tection scheme. The Standard Quantum Level (SQL) of 
the shot noise is obtained by subtracting their photocur¬ 
rents. The overall detection efficiencies accounting for 
photodetector efficiencies and losses in the optical paths 
are 87% for the twin beams and 65% for the pump beam. 
Electronic noise is subtracted according to independent 
calibration of its Gaussian distribution. 

The photocurrent is spectrally analyzed by a demodu¬ 
lating module, where the temporal signal is electronically 
mixed with a sinusoidal reference (eLO) at 12 = 21 MHz, 
in this manner defining the frequency of spectral field 
modes of interest. The result of electronic mixing is 
filtered by a 600 kHz bandwidth low-pass filter and 
recorded by a computer with an analog-to-digital board 
(NI PCI-6110). For each combination of field quadratures 
on a single beam, 1,000 spectral quantum measurements 
are realized in order to obtain photocurrent moments. 
In total, 450 different directions in the two-mode phase 
space of each beam are probed with RD, providing a to¬ 
tal of 450,000 spectral quantum measurements per beam 
in 750 ms acquisition time. 

FIG. 2. Experimental setup. Nd:YAG Diabolo: Pump laser; 
OPO: Optical parametric oscillator; PBS: Polarizing beam 
splitter; FR: Faraday rotator; KTP: Nonlinear crystal. 

Experimental results - We now examine the quan¬ 
tum statistics of photocurrent fluctuations of the tri¬ 
partite system composed of the pump, signal, and idler 
light beams produced by the above-threshold OPO. Our 
aim is to provide a detailed study of the Gaussian or 
non-Gaussian character of the resulting quantum states, 
which involve in principle 6 spectral modes (a pair for 
each beam) [l3j]. Different quantum state can be pro¬ 
duced by controlling the input pump power P = Pq/P th 
(measured relative to the threshold power). For each 
of them, we probe with RD the quadrature probability 
distribution in a given direction in the phase space of 
quadrature observables. 

We start by analyzing the quantum states of each in- 
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dividual beam (2 spectral modes) for different values of 
pump power P. For each quantum state (i.e. each value 
of pump power), we perform a thorough investigation 
of a single quadrature direction in phase space. Focus¬ 
ing on a single quadrature observable allows us to gather 
larger statistical samples and analyze the Gaussian char¬ 
acter of the photocurrent up to the 14 th -order moment 
(n = 7). The chosen quadrature corresponds to the semi- 
classical amplitude quadrature and is measured with the 
analysis resonators far off resonance (A 1). We col¬ 
lect N = 280,000 consecutive quadrature quantum mea¬ 
surements for each quantum state. One such data set, 
for P = 1.72, is presented in Fig. [3] From top to bot¬ 
tom on the left column, idler, pump, and signal spectral 
photocurrent fluctuations are shown. The corresponding 
histograms are presented on the right column of Fig. [3] 
They show clear visual agreement with the normal dis¬ 
tribution. 
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FIG. 3. (a) Raw data of quadrature amplitude fluctuations 
obtained with the analysis resonator far off resonance, (b) 
Histograms of raw data with the normal density distribution 
superimposed. 

Quantitative analysis of the Gaussian character of 
those signals is first performed by considering the 3 rd - 
and 4 th -order moments, shown in Fig. [I] We define the 
third order moment dj = a j 3 ^/s d , related to the skewness 
of the probability distributions, and fourth-order moment 
kj = erj 4 ^/s d , akin to the kurtosis. Indices j = 0,1, 2 in¬ 
dicate pump, signal, and idler beams, respectively. For 
each value of pump power, moments are extracted from 
the complete set of 280,000 data points. For Gaussian 
statistics, one expects dj = 0 and kj = 3, as indicated by 
the black dashed lines on the figures. 

According to the results in Fig. all measured fourth- 
order moments are compatible with the Gaussian distri¬ 
bution, meaning that the photocurrent signal of any sin¬ 
gle beam produced by the above-threshold OPO would 
fulfill Eq. (1121) within the experimental uncertainty. In 
the case of dj , the null result does not provide informa¬ 
tion about the quantum state, since it is a consequence of 
the phase mixing regime [Eq. (TTU1) ] . For both quantities 
the fluctuations observed in the datasets are compati¬ 
ble with the statistical uncertainty of each data point. 
The number N of quantum measurements per data trace 
alone would allows us to establish the Gaussian character 
of the photocurrent within 1 / y/N ss 0.2% of the moment 
values. However, at these levels of confidence, we could 


notice by the observation of higher than expected dis¬ 
crepancy between different experimental data points the 
presence of a systematic error affecting the mean value 
of the photocurrent. This error source, probably caused 
by the photodetection electronics and still under inves¬ 
tigation, has been included on the overall uncertainty of 
our data. The overall values and uncertainty for pump, 
signal, and idler fourth-order moments considering the 
whole data set as one are respectively k 0 = 2.9987(17), 
k\ = 2.9951(26), and fc 2 = 2.9931(27). The third-order 
moments are do = —0.0006(9), d\ = —0.0041(21), and 
d 2 = -0.0033(23). 



FIG. 4. Statistical analysis of 3 rd - (top) and 4 th -order (bot¬ 
tom) moments of single-beam photocurrent measurements rel¬ 
ative to the distribution standard deviation. Dashed lines 
represent the values expected for Gaussian statistics. Green 
circle: pump; Blue up-pointing triangle: signal; Red down¬ 
pointing triangle: idler. 

Moments of higher order are now probed in the same 
way, by keeping the same choice of quadrature measure¬ 
ment (A> 1). Results are presented in Fig. [5] We con¬ 
sider the ratios 7 'f 2n f = o-j 2 "^/s 271 , for n = 2,3,4, 5,6, 7, 
where Sj is the standard deviation of the quadrature 
probability distribution for beam j. Values expected for 
a Gaussian distribution are indicated by the dashed lines. 
All moments show perfect compatibility with the Gaus¬ 
sian distribution within the experimental uncertainty. 

Next we investigate whether different directions in 
phase space could show deviations from the Gaussian 
statistics, by varying resonator detunings. We consider 
the quantum state produced by a fixed pump power P 
= 1.66(5) and perform quadrature measurements with 
RD (similar results are obtained for different values of 
P). Fig. [6] depicts RD the spectral photocurrent of the 
three beams measured as functions of resonator detun¬ 
ings. Formally, each value of A corresponds to a direc¬ 
tion of observation in the two-nrode phase space. In¬ 
set (a) shows the measured quantum fluctuations for the 
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FIG. 5. Statistical analysis of higher-order moments of single¬ 
beam photocurrent measurements relative to the distribution 
standard deviation. Dashed lines represent the values ex¬ 
pected for Gaussian statistics. Green circle: pump; Blue up- 
pointing triangle: signal; Red down-pointing triangle: idler. 


pump beam as functions of A and normalized by its av¬ 
erage intensity (raw data), corresponding to single mea¬ 
surements of the observable Ig( A). Inset (b) presents the 
quantum noise (variance of photocurrent fluctuations) for 
the three single beams, relative to the SQL. Each point 
in the graph is a realization of Eq. m, averaged over 
1,000 single quantum measurements. While the pump 
beam presents nearly shot-noise limited quantum noise 
for all values of A, signal and idler beams show excess 
noise in some directions in phase space (in the semi- 
classical picture, the beams would be described as pos¬ 
sessing phase noise), determined by a theoretical curve 
used to obtain the optimum values of mixed operator 
moments. Fig. [G] (c) presents the results of kj for differ¬ 
ent quadratures of each single-beam, obtained by com¬ 
plete scans in A. The statistical uncertainty in this case 
is around 1/vlOOO « 3%, and no significant deviation 
from the Gaussian statistics is found. Inset (d) displays 
once more a null value of dj, a consequence of the phase 
mixing regime. 

Finally, the analysis of two-beam correlations, in which 
four spectral modes are involved, is realized by changing 
modal basis to the sum or subtraction of single beam 
spaces. We consider the operators ) = (/L) -j- 

ijjj ^)/-\/2, j 7 ^ j', constructed by the coherent combi¬ 
nation of single-beam observables. In this manner, two- 
beam probability distributions are probed with the tools 
employed above for single-beam quadratures. In fact, 
the presence of non-Gaussian correlations among beams 
would appear in this picture as non-Gaussian probability 
distributions of single quadrature observables related to 
two beams, by a change of modal basis. Fig. 0 summa¬ 
rizes the experimental results. No deviations from the 
Gaussian statistics can be observed for either the sum or 
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FIG. 6. Statistical analysis as a function of resonator detun¬ 
ing. (a) Quantum fluctuations of the pump beam normalized 
by the average intensity (raw data), (b) Spectral quantum 
noise produced by pump (green curve, bottom), signal (blue 
curve, middle) and idler beams (red curve, top), (c) Fourth- 
order moments normalized by the standard deviation, (d) 
Third-order moment normalized by the standard deviation. 
Dashed lines represent the exact value expected from Gaus¬ 
sian statistics. In (c) and (d), symbols follow the same code 
of Fig. [5] 


the subtraction of single-beam photocurrents in the data. 

We have also applied the Shapiro-Wilk normality test, 
or W test, to further confirm the Gaussian character of 
the photocurrent signals. The test determines the de¬ 
gree of reliability by which the data sample would be in- 
congruent with Gaussian statistics (the null-hypothesis). 
Shapiro and Wilk proved the power of their test show¬ 
ing its capability of detecting non-normality in small sets 
for a wide variety of statistical distributions, including 
those with Gaussian kurtosis values 22], The test can 
be extended to sets containing up to thousands of sam¬ 
ples. We have performed numerical investigations and 
found no advantage in the W test applied to our data, for 
which approximately 1% accuracy is reached. Due to our 
large sets of data (hundreds of thousands), the analysis 
of higher order moments presented above reaches better 
precision. 


IV. CONCLUSION 

When employing the spectral analysis of the photocur¬ 
rent to retrieve properties of the quantum state of light, it 
is important to realize that the Gaussian or non-Gaussian 
character of the observed photocurrent noise does not 
establish by itself the type of statistics followed by the 
quantum state, since the photocurrent formally results 
from an incoherent mixture of two independent quan¬ 
tum measurements. Pure photocurrent operators, better 
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FIG. 7. Statistical analysis of higher order moments of two- 
beam photocurrent measurements normalized by the standard 
deviation of the respective probability distribution. Dashed 
lines represent the exact value expected from Gaussian statis¬ 
tics. Sum (top) and subtraction (bottom) of single-beam 
photocurrents are shown. Green circle: pump-idler combi¬ 
nation; Blue up-pointing triangle: pump-signal combination; 
Red down-pointing triangle: signal-idler combination. 


understood in the modal basis of symmetric (5) and anti¬ 
symmetric (.A) combinations of spectral sideband modes, 
can only be retrieved if both the optical and electronic lo¬ 
cal oscillators are phase coherent to one another during 
the whole process of quantum state reconstruction, in an 
implementation of a doubly phase-coherent detection. 

Thus the observation of Gaussian statistics of the spec¬ 
tral photocurrent does not provide an unambiguous ac¬ 
count of the quantum state statistics with the usual 
incoherent detection. It either strongly constrains the 
types of non-Gaussian quantum states capable of pro¬ 
ducing the Gaussian photocurrent or provides strong ev¬ 
idence for Gaussian and symmetric single-mode quan¬ 
tum states of S and A modes. A similar scenario holds 
in case non-Gaussian photocurrent fluctuations are ob¬ 
served. While it is true that non-Gaussian quantum 
states could generate the observed statistics, an inco¬ 


herent mixture of Gaussian processes may also produce 
non-Gaussian quantum noise. 

The use of the photocurrent statistics to reconstruct 
the quantum state in general requires a priori knowl¬ 
edge. In most experiments in quantum optics, at least 
one optical field with ‘classical’ characteristics partici¬ 
pates in the quantum dynamics, in which case a linearized 
interaction will favor the interpretation of data as stem¬ 
ming from Gaussian and symmetric quantum states in 
the S/A modal basis (if one of those characteristics is 
taken as assumption, the other follows as consequence). 
In this scenario, the Gaussian character of the photocur¬ 
rent implies that the phase-locked photocurrent compo¬ 
nents (i.e. the ones which would be measured with dou¬ 
bly phase-coherent detection) are themselves stationary. 
Thus, even though phase mixing does not allow direct 
verification of this property, the Gaussian character of 
the mixed signal indirectly establishes it, as long as the 
quantum state is assumed to be Gaussian (or, equiva¬ 
lently, the possible non-Gaussian quantum states that 
would lead to the same Gaussian photocurrent are dis¬ 
missed as implausible). 

The correct reconstruction of the quantum state heav¬ 
ily relies on this fact, since only then can the covariance 
matrix be written in the simple form of an effective single¬ 
mode Gaussian quantum state fl4| . Most experimental 
work in quantum optics employing spectral noise analy¬ 
sis assumes stationarity to hold without explicitly men¬ 
tioning it, a limitation that could be a problem for the 
realization of CV quantum information tasks on those 
systems. Another consequence of the fact that most ex¬ 
periments with bright light beams will favor the genera¬ 
tion of Gaussian states of the field lies in the interpreta¬ 
tion of non-Gaussian photocurrent, which should first be 
likely attributed to a mixture of two different Gaussian 
distributions, one pertaining to the pure quadrature dis¬ 
tribution of mode 5, and another to mode A. To further 
substantiate the non-Gaussian character of the quantum 
state, the mixing of two Gaussian processes must first be 
ruled out, for instance by showing that mixed Gaussian 
models fail to account for the data. However, the unques¬ 
tionable establishment of non-Gaussian properties of the 
spectral quantum state would require stronger evidence 
to be convincing, owing to its prized value as a resource 
in the realization of quantum information protocols. 

We apply these considerations to perform a thorough 
experimental investigation of the photocurrent noise from 
pump, signal, and idler beams produced by the above¬ 
threshold OPO. Our data establishes the photocurrent 
quantum noise of those beams as Gaussian within a broad 
range of investigated parameters. We analyze moments 
up to the 14 th order for individual beams and for cross¬ 
correlations. From a theoretical point of view, the OPO 
dynamics above the oscillation threshold is expected to 
be completely described by a linearized model (since only 
bright beams drive the interaction), from which Gaussian 
states of the field follow. The experimental observation 
of Gaussian photocurrent together with the well founded 
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theoretical assumption of Gaussian quantum states im¬ 
plies the symmetry of S and A local quantum states for 
each beam. Our results demonstrate the stationarity of 
the phase-coherent spectral photocurrent to justify the 
use of the semi-classical model of quantum noise and thus 
validate our previous analysis of multipartite entangle¬ 
ment in this physical system (hUrI ]. Even though the 
complete state of the three beams entails six modes for 
a given analysis frequency, it is legitimate to restrict the 
analysis to three ‘effective’ modes (choosing either set of 
the equivalent S or A modes), by a partial trace opera¬ 
tion. The missing correlations between S and A modes, a 
feature yet to be measured, can only increase the strength 
of quantum correlations in case all six spectral modes are 
considered. 

We conclude that in most experiments a priori knowl¬ 
edge is unavoidable to perform spectral quantum state 
reconstruction. A few exceptions are experiments which 
use both the LO and the eLO to generate the quan¬ 
tum state, as in Refs. 0II1- Such an improved sit¬ 
uation can be applied to all experiments in quantum 
optics dealing with spectral modes by phase-locking the 
optical and electronic local oscillators, e.g. by employ¬ 


ing sub-Hz linewidth lasers [23|, to realize doubly phase- 
coherent detection. This would ensure the phase coher¬ 
ence of all oscillators employed to extract the spectral 
quantum noise, making the quantum measurement asso¬ 
ciated with the spectral photocurrent a pure observable. 
We view such experimental improvement as essential to 
achieve assumption-free realization of quantum informa¬ 
tion protocols with CV spectral modes, particularly in 
the case of quantum protocols requiring measurements 
of pure observables to be used as feedback to control the 
quantum state. 
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